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If detected, dark matter particles, such as WIMPs and sterile neutrinos, will be the earliest relics we can 
study, the first coming from before Big Bang Nucleosynthesis, an epoch from which we have no data so far. Here 
we discuss how different pre-BBN cosmological models affect the properties of these particles, and how these 
properties could thus allow to determine the evolution of the Universe before BBN. 



1. Introduction 

Dark matter (DM) particle candidates, such 
as WIMPs, weakly interacting massive particles 
(also sterile neutrinos and axions) are produced 
before Big Bang Nucleosynthesis (BBN). This is 
an epoch from which we have no data. BBN is 
the earliest episode (200 s after the Bang, T ~ 0.8 
MeV) from which we have a trace, the abundance 
of light elements D, '^He and ^Li. In order for 
BBN and all the subsequent history of the Uni- 
verse to proceed as usual, it is enough that the 
earliest and highest temperature during the last 
radiation dominated period, the so called reheat- 
ing temperature Trh, is larger than 4 MeV [1]. 

The relic density and relic velocity distribu- 
tion before structure formation of WIMPs (and 
other DM candidates) depend on the character- 
istics of the Universe (expansion rate, composi- 
tion, etc.) before BBN. If these particles are ever 
found, they would be the first relics from that 
epoch that could be studied. Thus we will want 
to extract as much information about the Uni- 
verse at the moment these particles decoupled as 
we can. Will we be able to differentiate cosmo- 
logical parameters from particle physics param- 
eters that determine their relic density and relic 
velocity distribution? That is, will we be able to 
discriminate between different pre-BBN cosmolo- 
gies through studying the DM particles? To start 
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with, we need to know how large are the possible 
effects of different viable pre-BBN cosmologies on 
DM particle properties we could measure. 

2. Standard and non-standard pre-BBN 
cosmologies 

The argument showing that WIMPs are good 
DM candidates is old. The density per comoving 
volume of non relativistic particles in equilibrium 
in the early Universe decreases exponentially with 
decreasing temperature, due to the Bolzmann fac- 
tor, until the reactions which change the particle 
number become ineffective. At this point, when 
the annihilation rate becomes smaller than the 
Hubble expansion rate, the WIMP number per 
comoving volume becomes constant. This mo- 
ment of chemical decoupling or freeze-out hap- 
pens later, i.e. for smaller WIMP densities, for 
larger annihilation cross sections a. If there is 
no subsequent change of entropy in matter plus 
radiation, the present relic density is il.^'^h'^ ~ 
10^^" GeV~'^ / (av) , which for weak order a gives 
the right order of magnitude of the DM density 
(and a temperature T/.o. — TOj^/20 at freeze-out 
for a WIMP x of mass m^). 

This is a ballpark argument. When actually 
applied to particle models, the requirement that 
the WIMP candidate of the model must have 
the measured DM density is very constraining. 
In many supersymmetric models, in which the 
WIMP candidate is usually a neutralino, this 
"DM constraint" is very effective in restricting 
the parameter space of models. With the stan- 
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Figure 1. Standard neutralino relic density 
vs mass for 1700 different MSSMs (colo indi- 
cate composition: green for bino-like, pink for 
higgsino-like, brown for wino-like) [6]. 

dard pre-BBN cosmological assumptions, after 
the constraints imposed by LEP II, having a 
neutralino with the DM density requires special 
tuned relation among parameters [5] . In Minimal 
Supersymmetric Models (MSSMs) that could be 
found at the LHC, the "DM constraint" requires 
a "well tempered neutralino" 5 , fine tuned to be 
at the boundary between a pure bino and a pure 
higssino or a pure bino and a pure wino. In fact, 
as shown in Fig.[Tl with the standard cosmological 
assumptions bino-like neutralinos tend to be over- 
dense (except when they annihilate resonantly or 
coannihilate with an almost degenerate slepton) 
while higgsino-like and wino-like are underdense 
(or with masses close to ITeV or 2 TeV respec- 
tively, beyond the reach of the LHC). Fig.[T]shows 
the relic density as function of the mass of neu- 
tralinos of different composition in 1,700 different 
MSSMs characterized by nine parameters defined 
at the electroweak scale [6]. In minimal super- 
gravity models (mSUGRA) for instance, the neu- 
tralino is typically bino-like thus its relic density 
tends to be larger than observed. The "DM con- 
straint" is found to be satisfied only along four 
very narrow regions in the fermionic and scalar 
mass parameter space mi/2, (see e.g. Ref. 1). 
Most of the "benchmark points", special models 
chosen to study in detail in preparation for the 



LHC and the next possible collider (such as A' to 
L', Snowmass Points and Slopes or SPS la', lb, 
2, 3 ,4, 5, Liner Collider Cosmo points or LCC 
1,2,3,4) lie on those very narrow bands 3 . Neu- 
tralinos are underabundant (account for a frac- 
tion of the DM) also in narrow regions adjacent 
to these just mentioned, but in most of the pa- 
rameter space neutralinos are overabundant and 
the corresponding models are many times said to 
be rejected by cosmology. Is it correct to reject all 
these supersymmetric models? The answer is no. 
The issue is that the position in parameter space 
of the narrow bands just mentioned depends not 
only on the particle model to be tested in collider 
experiments, but on the assumptions made about 
the history of the Universe before BBN. 

The standard computation of relic densities re- 
lies on assuming that radiation domination began 
before the main epoch of production of the relics 
and that the entropy of matter and radiation has 
been conserved during and after this epoch. In 
the case of WIMPs it is also assumed that they 
are produced thermally, i.e. via interactions with 
the particles in the plasma. With these assump- 
tions chemical decoupling happens at T^.o. — 
TO^/20 and kinetic decoupling, the moment after 
which WIMPs do not exchange momentum effi- 
ciently with the cosmic radiation fluid, happens 
later, at T^a. ^ 10 MeV-10 GeV [4]. Thus all 
WIMPs with > 100 MeV decouple at tem- 
peratures higher than 4 MeV, when the content 
and expansion history of the Universe may differ 
from the standard assumptions. 

Non-standard pre-BBN cosmological models 
are more complicated than standard ones and al- 
though several aspects of them have been stud- 
ied, a comprehensive scenario in each of them is 
yet missing. Usually non-standard cosmological 
scenarios contain additional parameters that can 
be adjusted to modify the WIMP relic density. 
However these are due to physics that does not 
manifest itself in accelerator or DM detection ex- 
periments. 

In non-standard cosmological models, the 
WIMP relic abundance may be higher or lower 
than the standard abundance. The density may 
be increased by creating WIMPs from decays of 
particles or extended objects (non-thermal pro- 
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Figure 2. H(T) for several pre-BBN cosmological 
models: radiation-dominated (RD), LTR, kina- 
tion (K) and scalar tensor with H increase (STl) 
and with H decrease (ST2). Fig. from Ref.[ig. 

duction) or by increasing the expansion rate of 
the Universe at the time of freeze-out. The den- 
sity may be decreased by reducing the expansion 
rate of the Universe at freeze-out, by reducing 
the rate of thermal production (through a low 
TfiH < Tf,o.) or by producing radiation after 
freeze-out (entropy dilution) . Any extra contribu- 
tion to the energy density of the Universe would 
increase the Hubble expansion rate H and lead to 
larger relic densities (since the decreasing inter- 
action rate becomes smaller than H earlier, when 
densities are larger). This can happen in Brans- 
Dicke- Jordan [7] model, models with anisotropic 
expansion [ 51715] , scalar-tensor (ST) |10lllll2ll3j 
or kination ^14 9^ models and other models [T5] . 
These models alter the thermal evolution of the 
Universe without an extra entropy production. 
Kination |14j is a period in which the kinetic 
energy ~ 0^/2 of a scalar field (j) (maybe a 
quintessence field) dominates. The relic density 
of WIMPs which decouple during kination may 
be increased by a factor ^/r70lO^(m^/lOOGeV), 
where rj^ = p^/ p-y < 1 is the ratio of the (j) ^-nd 
photon densities at T ~ 1 MeV. ST models of 
gravity have a scalar field coupled only through 
the metric tensor to the matter fields. The ex- 
pansion of the Universe drives the scalar field to- 
wards a state where the theory is indistinguish- 



able from General Relativity at T < T^, before 
BBN. Fig.dlshows H{T) for several pre-BBN cos- 
mological models. WIMPs freeze-out before T^, 
when H is larger, and although after there 
is a "reannihilation phase" , the net effect is that 
WIMPs densities are increased by factors of 10 to 
10^ PJJ. In some ST models H may be decreased, 
leading to relic densities smaller than standard by 
as much as a factor 0.1 [13 . The factors just men- 
tioned, 10^ to 0.1, are not large enough to bring 
any MSSM neutralino, for example, to have the 
DM density. 

Not only the value of H but the dependence 
of the temperature T on the scale factor of the 
Universe a is different than the usual a ~ T~^, 
if entropy in matter and radiation is produced. 
This is the case if a scalar field oscillating 
around its true minimum while decaying is the 
dominant component of the Universe just before 
BBN, for which ~ ~ a'^/"^. Models of this 
type include some with moduli fields, either the 
Polonyi field 117)18] or others [19], or an Affleck- 
Dine field and Q-ball decay [20], and thermal in- 
flation [21]. Moduli fields correspond to flat di- 
rections in the supersymmetric potential, which 
are lifted by the same mechanisms that give mass 
to the supersymmetric particles of the order of a 
few to lO's of TeV, and they usually have interac- 
tions of gravitational strength. The decays of the 
(j) field finally reheat the Universe to a low reheat- 
ing temperature Trh =i 10 MeV(m0/lOO TeV)^/^ 
(m^ is the (f> mass), which could be not much 
larger than 4 MeV. In these low temperature re- 
heating (LTR) models there can be direct pro- 
duction of DM relics in the decay </> which in- 
creases the relic density, and there is entropy 
generation, through the decay of into radia- 
tion, which suppresses the relic abundance. The 
change in relic density in LTR models is larger 
than in other non standard models. Both thermal 
and non thermal production mechanisms in LTR 
have been dis cussed [2:jl7ll7ll8l23ll9l:j4l25l26l271 
12812916130131] . mostly in supersymmetric mod- 
els where the WIMP is the neutralino. The 
WIMP relic density depends only on two addi- 
tional parameters (determined by the model at 
high energy scales) besides the usual ones: Trh 
and b/m^, the ratio of &, the net number of 
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Figure 3. LRT neutralino relic density O/i^ vs 
mass for different values of Trh and rj for the 
same 1700 MSSMs as in Fig. [T] 0. 

WIMPs produced on average per (j) decay and the 
(j) mass (29j . 6 is a highly model dependent pa- 
rameter |19l25l29j . Fig. [3] shows the LTR relic 
density of neutralinos in the same 1,700 MSSMs 
(the same whose standard relic density is shown 
in Fig. [T]) , each shown as one point in each panel. 
Fig. [3] shows that all points can be brought to 
cross the right DM density cyan line with suited 
combinations (in general not unique) of Tjjh and 
T] = b lOOTeV/m^. This means that neutrali- 
nos can have the DM density in (almost) all su- 
persymmetric models, provided the right values 
of TfjH and rj can be obtained (the exception 
being severely overabundant or underabundant 
very light neutralinos, rarely encountered in su- 
persymmetric models [51). This has important 
implications not only for colliders but for direct 
and indirect DM searches as well. For example, 
the region of viable supersymmetric models to be 
searched for in direct detection experiments ex- 
tends from well under 1 GeV to 10 TeV in neu- 
tralino mass [32] . 

Not only the relic density of WIMPs but 
their characteristic speed before structure for- 
mation in the Universe can differ in standard 
and non-standard pre-BBN cosmological mod- 
els. If kinetic decoupling happens during the re- 
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Figure 4. Mass within the kinetic decoupling 
horizon and free-streaming volumes {Md and 
Mfs, respectively) of Ultra Cold WIMPs with re- 
spect to their standard values as function of the 
standard kinetic decoupling temperature for dif- 
ferent reheating temperatures Tj^h [33| . 

heating phase of LTR models, WIMPs can have 
much smaller characteristic speeds, i.e. be much 
"colder" 33J , with free-streaming lengths several 
orders of magnitude smaller than in the standard 
scenario (see Fig. U). Much smaller DM struc- 
tures could thus be formed, a fraction of which 
may persist as clumps within galactic haloes and 
be detected in indirect DM searches. The sig- 
nature would be a much larger boost factor of 
the annihilation signal than expected in standard 
cosmologies for a particular WIMP candidate. 
WIMPs may instead be much "hotter" than in 
standard cosmologies too, they may even be warm 
DM instead of cold, which would leave an imprint 
on the large scale structure spectrum [51] . 

Let us add that also finding "visible" sterile 
neutrinos, i.e. those that could be found soon 
in neutrino experiments, would be a signature of 
a non standard pre-BBN cosmology. If they are 
produced through oscillations, they are produced 
mostly at T ~ 13 MeV(ms/l eV)^/^, which is >4 
MeV if their mass is > 10"'^ eV . In order to 
be found in experiments, these sterile neutrinos 
would necessarily have mixings with active ones 
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large enough to be overabundant, and thus be 
rejected, in standard cosmologies. In LTR mod- 
els the relic abundance of visible sterile neutrinos 
could be reduced enough for them to be cosmo- 
logically acceptable [36] , 

3. Colliders as DM and pre-BBN cosmol- 
ogy probes 

In most scenarios one can think of the LHC 
should find at least a hint of new physics, and 
whatever it finds will lead to a set of possible DM 
candidates and reject others. In our most am- 
bitious scenarios, the LHC will be able to deter- 
mine a range for the annihilation cross section of 
the DM candidate, thus a range for the standard 
relic density Q^'^. For example, the standard relic 
density in the LCC2 model could be determined 
within 40%, and the error could be reduced to 
14% with the ILC-500 GeV and to 8% with the 
LHC- 1000 GeV [3S]. 

Assume that the LHC finds a DM candidate 
and determines its standard relic density to be 
much larger than the DM density. Then, either 
the particle is unstable (e.g. in SUSY models it is 
the next-to-lightest particle, the NLSP, and the 
SUSY spectrum should tell if this is true) or it 
does constitute the DM (found in DM searches) 
and the pre-BBN cosmology is non standard. 

The LHC may instead determine that the stan- 
dard relic density of the DM candidate is much 
smaller than the DM density. In this case there 
are at least two possibilities: either this is just a 
part of the DM and the DM has other components 
(which are to be found through DM detection ex- 
periments), or this particle constitutes the bulk 
of the DM (and this will need to be determined 
through DM searches) and the pre-BBN cosmol- 
ogy is non standard. The separation between 
these two possibilities will be difficult, in part be- 
cause the interaction rate in direct DM detection 
experiments may be large even for very under- 
abundant DM halo particles [37]. Indirect DM 
searches may be very important in this case |38| . 

If ^l^'^ of a DM candidate is found to be com- 
patible, within the error in the determination, 
with the DM density, there are still at least two 
possibilities. One is that the DM candidate found 



is actually not the DM (namely it is not found 
in DM searches even if the interaction cross sec- 
tion is large enough). In this case either the DM 
candidate found is unstable and decays outside 
the detectors, or the pre-BBN cosmology is non- 
standard and the particle has a very small relic 
density. The second possibility is that we have 
truly found the main component of the DM. Even 
in this case we still would want to get bounds on 
the possible departure of the pre-BBN cosmology 
from standard, given that this particle is the only 
relic we have from this early epoch in the history 
of the Universe. This, in fact, may be a good 
argument (one among many) to build the next 
collider after the LHC. 

In summary, so far we have explored how well 
the LHC and ILC could determine the relic den- 
sity largely assuming that the pre-BBN evolution 
of the Universe is standard. In reality, with the 
LHC and ILC we are trying to measure a com- 
bination of both, the DM relic density and key 
parameters of the pre-BBN cosmology which are 
necessarily tied up. To disentangle both types of 
parameters we need to combine accelerator mea- 
surements with direct and indirect DM searches, 
which would tell us the actual halo density of the 
DM particle candidate. 

The best possible calculation of the WIMP relic 
density (and relic velocity distributions) assuming 
a standard cosmology are necessary to eventually 
discriminate between standard and non standard 
pre-BBN cosmologies, but not to reject elemen- 
tary particle models at this time (no "DM con- 
straint" should be used to reject particle models). 

Few attempts have been made so far to figure 
out how to discriminate particle physics parame- 
ters and pre-BBN cosmological parameters at col- 
liders. For example, Ref. [39] studied benchmark 
points LCCl', 2', 3' and 4', slightly shifted in pa- 
rameter space with respect to the unprimed LCC 
points, so that neutralinos are underabundant in 
the standard cosmology. Assuming neutralinos in 
these models decoupled while the Universe was in 
a kinetion period, a simultaneous search could be 
made at the LHC and ILC for the relic density 
and for the kinetion parameter rj^. Ref |39j finds 
that in some of these models the LHC can distin- 
guish that the relic density is the DM density for 
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T]^ non-zero, but in others the LHC finds that the 
models are compatible with the standard cosmol- 
ogy. In all cases the ILC would do a good job of 
determining both parameters. 

Much more work is needed to try to figure out 
how to discriminate particle physics parameters 
and pre-BBN cosmological parameters using the 
data of colliders and other experiments on DM 
candidates. 
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